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THE CALCULATION OF SEPARATED FLOW AT HELICOPTER BODIES 


G. Polz 

Messerschmidt-Bolkow-Blohm GmbH, Munich, W. Germany 

1. InCroduction 

In the design of helicopter fuselages, designs frequently 
result — due to the end-use of the aircraft — which are not optimal 
from an aerodynamic standpoint (fig. 1). For these fuselage shapes, 
flow separations frequently occur (fig. 2) which are noticed pri- 
marily by the resistance, and thus indirectly by the economics, 
and also by a deterioration in stability properties and vibrations 
with similar effects. Due to the observed trend in modern heli- 
copters toward higher flight speeds, this problem — which was of 
subordinate importance to earlier developments — has recently moved 
to the foreground. The problems connected with this provided the 
impetus for various investigations on flow separation from heli- 
copter fuselages, which were partly of an experimental nature to 
solve real problems [1,2], and partly of a theoretical nature for 
basic research on physical relationships [3-12]. 

2. Flow Sepai'ation and Wake Formation 

Among the various forms of flow separation (fig. 3), for the 
typically compact fuselage shapes of helicopter, wake separation 
occurs, whereas eddy separation occurs only at extreme angles of 
attack. The reason for the wake separation is a relatively pre- 
cipitous decrease in cross-section in the fuselage region and 
the retarded flow connected with it. The continuation of the 
boundary layer forms the wake region which contains a backflow 
region near the wall. 

Since the flow material coming from the boundary layer is 
rotational or frictional, the wake contains no potential flow, 
which impedes a computational treatment of such flows. The rotation 
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present there corresponds in its effect, to a vortex distribution 
whose induction affects the entire flow field around the fuselage. 
Thus the description of wake influences is of prime importance in 
the determination of this induction effect. 


The origination of wake circulation can be illustrated on the 
example of a rotation-symmetric body with incident flow (fig. 4): 
The entire wall boundary-layer represents a shear-flow region — 
since the velocity jump from u = 0 at the surface to the outside 
speed u = U must be overcome in it — for which the Euler equation 
for rotational flows can be used: 
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( 1 ) 


In a local coordinate system, the lateral component v and all 
derivatives of y disappear due to rotational symmetry; in addition 
3w/8x can be neglected. Thus it follows that the rotation 

u = ^ (0, 3u/3z, 0) (2) 


depends only on the velocity gradient. The circulation 

, 3r 6 '5 a 5 

= ^ = 2 / m • dz = / ^ dz = / du = U 

y 3 s 0 y 0 0 ( 3 ) 

present in the boundary layer per length unit depends only on the 

outside velocity U, but not on the velocity profile or on the thick- 
ness of the boundary layer. Since the circulation in a boundary 
layer element of length ds 

dr^ = j^V»dl = d({) = U*ds (4) 

is equal to the change in potential d({), this means that each vortex 
element is bounded by two equipotential lines. From the arrangement 
of potential lines we can see that the entire circulation present 
in the boundary layer can be described by a system of closed eddy 
lines oriented in the circumferential direction. 
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The vortex intensity streaming into the wake with the boundary- 
layer material at the separation line 
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(5) 


represents the sole source for the "circumferential" circulation 
present in the wake. The backflow present in this region likewise 
leads to the formation of a boundary layer whose rotation is 
opposite that of the main boundary layer. 


In general however, non-rotation-symmetric incident-flow 
states are present, connected with a cross-force (e.g. lift) acting 
on the body. This cross-force causes the appearance of natural 
vortexes whose free ends are parallel to the incident flow direc- 
tion. In contrast to the typical vortex separation, these vortex 
filaments are present in the compact helicopter fuselage as uni- 
formly distributed "longitudinal" circulation within the wake 
(fig. 2). The wake contains both a circulation oriented in the 
circumferential and in the longitudinal direction (fig. 5). For 
the generation of longitudinal circulation, the direction of flow 
lines at the fuselage surface directly in front of the separation 
line is probably decisive. 


3. Previous Methods for Computing the Separated Flow from 

Helicopter Fuselages 

Up until a few years ago, no computing methods were known for 
the calculation of the separated flow from general fuselage-shapes 
like helicopters. But recently, the development of methods was 
taken in hand specifically for helicopter fuselages, including 
efforts in the USA [3,4,5]. 

\ 

3.1 Simulation of the Wake by Flow from the Fuselage 

An initial test was performed by Woodward [3] of the fuselage 
surface lying in the wake region by using an additional source 
distribution (for the fuselage itself the panel method [14] based 
on a source distribution of the surface was already used). A 
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similar method was developed by Jacob [13] for computing the 
separated flow from profiles. This method does permit a sufficiently 
accurate determination of the pressure distribution on the fuse- 
lage surface, but the pulse depression and thus the flow condi- 
tions in the far wake thus cannot be reproduced with it. A sat- 
isfactory agreement with measured results could only be achieved 
(fig. 6) when only a part of the fuselage surface lying in the 
separation region was occupied by additional sources. The source 
strength was determined for each of the affected panels by the 
equation: 


V. 


n = V 


n 


( 6 ) 


i.e. the normal component of the velocity induced by the source 
was set equal to the normal component of the incident-flow velocity, 


3.2 Wake Simulation with Vortex Distribution 

Another method which permits a more realistic simulation of 
the wake was developed by Dvorak et al. [4,5]. The fuselage is 
again represented by panels and a vortex surface lying parallel 
to the incident-flow direction is first set at the panel edge 
lying closest to the separation line (fig. 7). The circulation 
of this vortex surface is determined by a control point on that 
surface. The length L of the vortex element* is selected so that 
the velocity along the vortex surface remains constant in a certain 
range, which requires an iterative computation of the final pressure 
distribution. Characteristic of the results is an excessively 
stressed pressure increase in the region in front of the point 
of separation. Otlierwise, a relatively accurate reproduction of the 
measured pressure distribution is possible. 

4. Development of a Computational Model with Uniform Vortex 

Distribution in the Wake 

Within the frame of a ZTL-contract , the development of an 
in-house computational method was begun at MBB for separated 
flows from helicopter fuselages, with the emphasis being placed on 
the simulation of flow states in the far wake. 
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4.1 Theoretical Model 


The singularity methods used in flow mechanics are all based 
on an integral formula of the form; 
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which cannot be solved for the general case. It can be solved 
only be breaking it down into individual elements, e.g. of the 
surface of the body. Since just as many boundary conditions as 
singularities are present, a closed solution to the integral 
equation is possible via a linear equation system. In the case 
of a separated flow, the left side of the equation is expanded by 
the influence of the circulation contained in the wake oj (fig. 8). 
This circulation cannot be determined by additional equations — in 
contrast to the case for an airplane wing — since neither its geo- 
metry nor its boundary conditions are known. But its calculation 
is possible by an iterative solution to the equation system by 
including in the boundary calculation, the induction of the wake 
circulation at each control point. Since this induction effect is 
not known from the beginning, it must be assumed in a suitable 
manner, for the first iteration step. Then with the first solu- 
tion for the singularity distribution, the geometry and circula- 
tion distribution of the wake can be determined — with additional 
appropriate assumptions. 

4.2 Structure of the Computational Method 

In accord with the presumptions of the theoretical model, 
the computational method consists of several program segments 
(fig. 9): The compression effect of the fuselage is computed with 
a panel method [14,15] which is based on a source distribution 
of the fuselage surface. By interpolation of the velocity distri- 
bution obtained with it, a number of flow lines is determined on 
the fuselage surface. For each of these flow lines, the separation 
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point and the boundary-layer thickness present there are deter- 
mined by a 2D-boundary-layer method [16]; in this case, the fuse- 
lage is replaced for each flow line by an equivalent rotation body. 
The line connecting the individual separation points represents 
the front edge of a trailing body which is set on the fuselage 
and covered by a defined circulation. The frame of this trailing 
body is formed by the free flow lines whose starting points are 
away from the fuselage surface by the thickness of the boundary 
layer and lie above the separation points of the flow lines. Since 
flow material is continually moved from the outside flow into the 
wake due to friction effects, the wake's cross-section continually 
increases compared to these flow lines. The circulation present 
in the trailing body again affects the singularity distribution 
of the fuselage. Since the influence of the v?ake is not known in 
advance, the calculation of the flow state requires an iteration 
process which uses a prismatic trailing body of constant cross- 
section for the ^first step. 

4.3 Geometry and Circulation Distribution of the Trailing Body 
Wake Geometry 

For the determination of the spread of the wake, the Prandtl 
mixing-path theorem can be used. The application of this theorem 
to the wake of different, blunt bodies [17] shows that the formulas 
based on it for the cross-sectional increase in a rotation-symmetric 
wake ( fig. 10) : 
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and the pertinent pulse indentation depth: 
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are valid for bodies of finite thickness only from a certain 
distance from the separation line, i.e. when the disturbance 
has already decayed due to the compression action of the body. 


There is no information available about the region near the 
body. Since the frictional wake-flow is continually enhanced by 
flow material from the outside flow only from the separation line, 
the wake cross-section must increase compared to the flow lines 
forming at the separation line toward the outside edge of the wake. 
For a numeric determination of cross-sectional development, the 
equation 



X - X. 


(1 + m) 


■A 


( 10 ) 


is selected; with m = 0.03 ... 0.05, it gives useable trailing-body 
shapes . 


Circumferential and Longitudinal Vortex Circulation at the Separation 
Line 

It seems useful to use, in accord with section 2, the velocity 
components present in the meridian-section to determine the cir- 
cumferential vortex, and the circumferential components to determine 
the longitudinal vortex. But unfortunately, the used boundary-layer 
computation method (applies precisely only for rotation bodies) does 
not always give the precise flow direction at the separation line, 
which is why a sufficiently accurate determination of longitudinal 
circulation is not possible in this case. On the other hand, when 
the zero-lift direction of the fuselage is known, the middleline 
of the trailing body which is sloped to the incident flow direction 
by ca. = (a - a©) /2 , can be determined. From the angle between 
the incident-flow direction and wake midline, a circumferential 
component can be determined (fig. 11) which leads to realistic 
values for the circumferential circulation. 
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Circulation Distribution in the Trailing Body 


A sufficiently accurate determination of the wake circulation 
is impossible due to the numerous unknowns (local flow direction, 
location and intensity of the individual vortices); thus one must 
rely mostly on empirical assumptions. At the separation line, the 
circumferential circulation is given by the difference between U = 0 
at the fuselage surface and the velocity prevailing at the edge 
of the boundary layer (fig. 12, position A). Under the assumptions 
of eq. (2), there results a parabolic circulation distribution w(z). 
Upstream from this (position B), a larger velocity difference 
is present due to the backflow, but this is again partly decreased 
by the wall boundary-layer of the backflow. For the circulation 
distribution this means that besides an intensive circulation in 
the outer layer of the wake, a weaker, opposite-rotating circulation 
is present near the wall. Behind the body there is generally a 
backflow region (C), v?hich is why Che velocity jump and thus the 
total circulation must be greater than at the separation line, 
assuming a nearly constant outside velocity U along the wake-mantle 
line. Outside the backflow region, the velocity difference and thus 
the circulation decreases gradually in accord with the Prandtl 
exchange law (D). A quantitative determination of the circumfer- 
ential circulation along the wake axis is impossible due to these 
considerations . 

Assuming that a similar velocity profile is present in the 
entire wake (with the exception of the fuselage region) as at the 
separation line — which is at least qualitatively correct — the same 
parabolic circulation distribution can be postulated everywhere. 

No quantitative statements are possible about the distribution 
of longitudinal circulation along the wake axis. But it can be 
asumed that it decreases gradually from the separation line due to 
frictional effects. For the radial distribution of longitudinal 
circulation, the same parabolic profile can be assumed as for the 
circumferential circulation. Otherwise, with respect to the circu- 
lation distribution of the trailing body, one must still rely on 
experiments . 
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Calculation of the Static Pressure 

The flow in the wake has a lower energy level than the outside 
flow, which is why the Bernoulli equation cannot be used here to 
compute the static pressure. This energy difference is computed 
from the velocity difference or from the circulation intensity Y 
of the eddy-layer which separates the considered field point from 
the outside flow; thus, the static pressure in the wake results as: 

Jr CO 

Numeric Calculation of the Wake Induction 

For this calculation, the spatially distributed circulation 
distribution of the irregularly shaped trailing body is replaced 
by suitable vortex elements. The geometric division of the trail- 
ing body defined by its mantle lines takes place through cut planes 
perpendicular to its longitudinal axis whose position in the body 
region is specified by the separation points of the individual 
flow lines, whereas behind the body, an equidistant cut-distribu- 
tion is selected. Each cut surface represents a polygon which is 
subdivided into individual triangles by radial rays to the vertex 
points (fig. 13). The ring-like circumferential vortex distribu- 
tion of the cut surfaces is replaced, for each triangle, depending 
on its distance to the particular field point, by one or more line 
vortices lying parallel to the outside edge of the triangle. The 
longitudinal vortex distribution of each triangular tube composing 
the trailing body, is viewed as constant. It too, is represented 
by one or more line vortices. To avoid singularity problems at 
field points lying close to or in the trailing body, additional 
cut planes are placed near the particular field point at a pre- 
cisely defined distance to the field point. For the same reason, 
the distribution of vortex elements in each triangle is specified 
in accord with the position to the field point. 

4.4 Results and Comparison with Measurements 

Figures 14 and 15 show the pressure distribution and the flow 
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field of a sphere at Re = 0.45 * 10°. With a circulation intensity 
constant at up to about one sphere-radius behind the separation 
line and then decreasing linearly, the best agreement was achieved 
here with the measured results. The pressure increase occurring 
in the calculation in the region of the separation point is 
attributable to the fact that the front edge of the trailing body 
acts like a step on the outside flow. An even better agreement 
with the measurements is obtained when the static pressure pre- 
vailing at the outside edge of the trailing body is used also for 
the sphere surface. The flow field (fig. 15) shows an at least 
qualitatively meaningful wake geometry with a backflow region. 

Figure 16 shows the location of the separation line at the 
BO 105 fuselage at different angles of attack by measuring and 
calculation with and without wake influence. Although no longi- 
tudinal circulation was taken into account here, the results (ob- 
tained after two iteration steps) generally agree with the measure- 
ments, at least for negative angles of attack. 

In figure 17 the fuselage flow lines and the trailing body 
including the pulse indentation are presented for the fuselage of 
the BK 117 at a = -5°. The deflection of the fuselage flow lines 
discernable in front of the separation line does not correspond 
to the profile obtained in the wind tunnel test (fig. 22), but 
the separation line is reproduced almost correctly. The trailing 
body consists of 15 flow lines; this figure has also proven useful 
in other calculations. There are no measurements of the velocity 
distribution in the wake, but the resistance computed from it 
corresponds rather accurately to the value measured in the wind 
tunnel. In the computation of the pressure distribution, there 
is a better agreement with measurements (fig. 18) than for the 
pure computation based on potential theory. The fluctuations in 
the pressure profile discernable in the BCL behind the separation 
line are caused by the local circulation distribution. 

Most recently, the method was used within the frame of an 
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investigation on optimurr! design of transport helicopter airframes 
[11] where the relatively detailed wake measurements were performed. 
Figure 19 shows the panel distribution for the wind tunnel model 
used in the calculation. In fig. 20, the profile of fuselage 
flow lines and the separation line from calculation and measurement 
are compared. Whereas the separation line is reproduced very well, 
the flow lines of the calculation are not sloped downward enough 
in the region of the separation line. With the longitudinal circu- 
lation defined in sec. 4.3, there results a velocity field in the 
wake cross-section lying directly behind the tail boom attachment, 
which comes very close to the measured one (fig. 21). 

4.5 Further Development of the Method 

Due to the unknown flow processes in the wake which are also 
almost always instationary, when applying the method an empirical 
adaptation of the individual parameters like wake spreading and 
circulation distribution, is needed. This is due primarily to the 
absence of suitable test results which would permit an adaptation 
of the computational model. But it is expected that the investi- 
gations of helicopter fuselages [11] begun last year in the DFVLR 
will permit a much better insight into the flow conditions pre- 
vailing in the wake. The results achieved to date represent 
considerable progress. 

In order to give an insight into the actual problems of flow 
separation in helicopters, two results from wind tunnel and (flight 
tests are provided: Figure 22 shows results from painting tests on 
the BK 117-model, both with and without spoiler on the fuselage 
rear. Due to the spoiler, the downward pressure present in the 
fast-flight state is reduced and the vortex pair connected with 
the downward pressure is moved downward which improved the flow 
conditions in the area of the tail unit. (In the meantime, the 
spoiler on the BK 117 has been replaced by another solution) . 

The computational description of the spoiler is connected with con- 
siderable problems, in part because of the complicated panelling. 
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Figure 23 comes from the flight testing of the BK 117 and 
clearly shows the various zones of separated flow by means of the 
wool threads. Besides the rear of fuselage where the separation 
is still intensified by the spoiler, above all the fuselage top 
behind the rotor shaft and parts of the tail surface are affected. 
The picture clearly shows that some development effort is still 
needed for a complete theoretical description of the separated 
flow at the helicopter airframe. 

5 . Summary 

The function-induced shapes and additionally occurring un- 
favorable incident-flow states are the principle causes for flow 
separations from helicopter fuselages. Due to these separation 
processes, both the resistance is increased and the effectiveness 
of. controls in the tail section is reduced. 

The processes of flow separation are explained, specifically 
with respect to the compact-shaped helicopter fuselage shapes, 
with emphasis on the circulation contained in the wake. In addi- 
tion, the existing methods where the wake simulation uses relatively 
simple means, are presented. Finally, a computational model is 
presented which operates with a uniform circulation distribution 
of the wake. The equation on which the method is based, the 
program run and the calculation of the influence of the wake are 
explained. The results include both pressure distributions and 
flow-line profiles and flow patterns and partly show a good agree- 
ment with measurements. On the basis of practical examples, the 
necessity of further development of such methods is demonstrated. 
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Figure 1: Typical Fuselage Shapes of Helicopters 



BK117 (a=ii»,3=o'’) 


Transport helicopter Model 
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Wake separation 


Vortex-layer Separation 


Figure 3; Wake and Vortex-Layer Separation 



Figure 4: Generation of Wake Circulation 

Key: 1-boundary Layer 2-separation line 3-wake 4-equipotential 
lines 
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Figure 5: 
Key: 1-wake 


Vortex Distribution of the Wake 
2-separation line 


z 



Fig. 6: Simulation of the Wake by Flow From the Fuselage Body 

Key: 1-without influence of wake 2-flow from A 3-flow from A and B 
4-measurement 
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^ PANEL MIT KONSTANTER QUELLSTARKE 



Fig. 7: Computational Model with Vortex Distribution of the 
Wake Outside Edge 

Key: 1-panel with constant source strength 2-panel with constant 
vortex strength 3-control point 4-measurement 5-without 
the influence of wake 6-iteration 7-fuselage, underside 
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Fig. 8: Numeric Solution to the Integral Equation for the Wake Flow 

Key: 1-integral equation 2-sum equation 3-non-linear equation 
system 
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Figure 9: Program Outline of the Computation Method for Separated 
Flow 

Key: 1-fuselage coordinates 2-geometric conversion 3-panel geometry 
4-sdlution of the potential equation 5-velocities and singular- 
ities on the body's surface 6-flow line computation 7-flow 
line profile on the body's surface 8-boundary layer calculation 
9-separation line 10-after-flow calculation 11-after-flow 
geometry, circulation intensities 



Figure 10: Geometry of the Trailing Body 
Key: 1-flow line 
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g. 11: DeteriTiLnat 
Vortex Cir 

y: 1-flow lines 
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g. 12: Velocity a 
y: 1-backflow reg 
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Figure 13: Calculation of Vortex Induction 


Fig. 
Key : 



14: Pressure Distribution on the-Sphere after Calculation and 
Measurement (Re = 0.45 • 10°) 

1-without separation (potential theory) 2-measurement 
3-calculation 4-pressure profile at the outside edge of the 
trailing body 




\\ 


Figure 15: Velocity Field of the Separated Flow on the Sphere 
(Re = 0.45 • 10^) 
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Figure 16: Separation Lines at the BO 105 Fuselage 

Key: 1-measurement 2-calculation with wake 3-calculation without wake 
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Fig. 17: BK 117 - Fuselage with Trailing Body 

Key: 1-separation line 2-flow lines 3-trailing body 
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Fig. 18: Pressure Distribution on the BK 117 - Wind Tunnel Model 
Key: 1-separation line 2-without wake 3-with wake 4-measurement 


5-calculation 
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Fig. 19: Panelling of the Transport Helicopter Model 



Transport Helicopter 
model 
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Fig. 20: Flow Profile at the Fuselage and Separation Line; Comparison 
of measurement and calculation 
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Figure 21: Velocity Field in Cross- 
the Transport Helicopter 
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